An approximate renormalization-group transformation for Hamiltonian systems with 

three degrees of freedom 

C. Chandre 1 , H.R. Jauslin 1 , G. Benfatto 2 , and A. Celletti 3 
1 Laboratoire de Physique- CNRS, Universite de Bourgogne, B.P. J^7 870, 
F-21078 Dijon, France 

2 Dipartimento di Matematica, Universitd di Roma "Tor Vergata", Via della Ricerca Scientifica, 

1-00133 Roma, Italy 

3 Dipartimento di Matematica Pura e Applicata, Universitd di L'Aquila, Via Vetoio, 

1-67100 L'Aquila, Italy 



<3\ 



in 



> 
oo 

o 
m 
o 

On 

6 

O 



We construct an approximate renormalization transfor- 
mation that combines Kolmogorov- Arnold- Moser (KAM) and 
renormalization-group techniques, to analyze instabilities in 
Hamiltonian systems with three degrees of freedom. This 
scheme is implemented both for isoenergetically nondegen- 
erate and for degenerate Hamiltonians. For the spiral mean 
frequency vector, we find numerically that the iterations of 
the transformation on nondegenerate Hamiltonians tend to 
degenerate ones on the critical surface. As a consequence, 
isoenergetically degenerate and nondegenerate Hamiltonians 
belong to the same universality class, and thus the corre- 
sponding critical invariant tori have the same type of scal- 
ing properties. We numerically investigate the structure of 
the attracting set on the critical surface and find that it is 
a strange nonchaotic attractor. We compute exponents that 
characterize its universality class. 



I. INTRODUCTION 

The breakup of invariant tori is one of the key mecha- 
nisms of the transition to chaos in Hamiltonian dynamics. 
For two dimensional systems and for quadratic irrational 
frequencies, it has been observed that, at the transition, 
a sequence of periodic orbits approaches geometrically a 
torus of the given frequency, with a nontrivial scaling be- 
havior This self-similarity has been described in 
terms of a nontrivial fixed point of a renormalization- 
group transformation in the case of the golden mean 
. The sequence of periodic orbits responsible for the 
breakup is generated by the continued fraction expan- 
sion of the frequency. For the extension to systems with 
three degrees of freedom (d.o.f.) involving three incom- 
mensurate frequencies, we lack a theory that generalizes 
the continued fractions. Numerically, three d.o.f. Hamil- 
tonian systems (or equivalently four dimensional volume- 
preserving maps) have been studied with an extension of 
Greene's criterion p~0|— |T^] and by reconstruction of invari- 
ant tori using conjugation theory j^M. The conclusion 
of these analyses was that there is no geometrical accu- 
mulation of periodic orbits around the critical torus, and 
thus absence of universality (at least for the specific fre- 
quency vectors they considered). 



In Ref. J17[ , an approximate renormalization-group 
scheme was described for a reduced family of isoener- 
getically degenerate Hamiltonians, which is an interme- 
diate case between two and three d.o.f., that in appropri- 
ate coordinates can be interpreted as one d.o.f. system 
driven by two periodic forces with incommensurate fre- 
quencies. This was also the class of models considered in 
Rcf. |Tl],[l2|,[l4| : They studied an intermediate case be- 
tween two-dimensional and four-dimensional volume pre- 
serving maps. In particular, invariant tori in these inter- 
mediate models act as barriers in phase space (limiting 
the diffusion of trajectories) as for two d.o.f. Hamiltonian 
systems. 

The conclusion of Ref. |l7j was that one can still ex- 
pect universal behavior in the breakup of invariant tori: 
The universality is associated to a hyperbolic nonpcriodic 
attractor of the renormalization flow. The idea is that 
all Hamiltonians attracted by renormalization to this set 
will display sequences of scaling factors that appear in 
a different order but with a universal statistical distri- 
bution. The dominant unstable Lyapunov characteristic 
exponent determines the approach to criticality of the 
universality class. 

By analyzing trajectories on the critical surface, we deter- 
mine the structure of the critical attractor. Our analysis 
indicates the existence of a strange nonchaotic attrac- 
tor whose correlation dimension seems to have a value 
around 1, but it is difficult to analyze it even numeri- 
cally. 

We take the definitions as formulated by Grebogi et 
al. : An attractor of a map is called strange if it is not 
a finite number of points, nor a piecewise differentiable 
set. An attractor is chaotic if typical orbits on it have a 
positive Lyapunov exponent. The strange attractor we 
obtain is of the same general type as the ones found in 
Ref. j[^,[l9), for quasiperiodically forced systems (of one 
and two dimensions). Some of their properties have been 
rigorously analyzed in Refs. |p0|-p2[. 
Chaotic attractors for renormalization maps have been 



conjectured and observed in statistical mechanics [£3 24 
and dynamical systems J25|-|2£| . In Ref. [^9],[50) , a strange 
chaotic attractor was found for renormalization of cir- 
cle maps, and in Ref. ]3l| l similar evidence was found 
for area-preserving maps, from the scaling analysis of 
periodic orbits. The origin of randomness in these lat- 
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ter studies is due to the randomness of the sequence of 
continued fraction approximants for an ensemble of the 
considered frequencies. In contrast, in the present three 
d.o.f. case, the rational approximants are a regular se- 
quence, obtained by iteration of a single unimodular ma- 
trix, which allows us to define a renormalization trans- 
formation with a fixed frequency vector. We observe that 
this difference leads to a qualitatively distinct structure 
of the attractors. Although the geometry of the attractor 
is singular (i.e., strange), it is not chaotic. 

In the present article, we extend the approximate 
renormalization-group transformation developed in Ref. 
jl7j to a more general family of Hamiltonians with three 
d.o.f. We find that the renormalization trajectories on 
the critical surface converge to the reduced family of 
Hamiltonians considered in Ref. 

Similar types of systems were first studied in Ref. ]32] ] 
with an approximate renormalization scheme that kept 
less terms than needed to detect an attractor. Instead, 
their transformation yields a rotation (a center) which is 
structurally unstable as mentioned in Ref. p2| . 

We define a renormalization transformation which acts 
on Hamiltonians with three d.o.f. written in actions A = 
(A 1 ,A 2 ,A 3 ) e JR 3 and angles ip = {(pi, (p 2 , ip 3 ) G T 3 (the 
3-climcnsional torus parametrized, e.g., by [0, 2ir} 3 ) 



H{A,<p)=H {A) + V(A,<p), 



(1.1) 



where Ho is the integrable part of the Hamiltonian. We 
are interested in the stability of the torus with frequency 
vector ljq. We suppose that this torus is located at A = 
for Hq, i.e., the linear part of Ho is equal to u>o ■ A. 
Kolmogorov-Arnold-Moser (KAM) theorems were proven 
for Hamiltonians (1.1) (with suitable restrictions on the 



perturbation) provided that Ho contains a twist in at 
least one direction in the actions p3j, i.e., the Hessian 
matrix d 2 Ho/dA 2 , with elements d^Io/dAidAj, is non- 
zero, and Wo satisfies a Diophantine condition. It shows 
the existence of the torus with frequency vector u)q for 
a sufficiently small and smooth perturbation V. The in- 
variant torus is a small deformation of the unperturbed 
one. To assume nonzero Hessian matrix, we restrict the 



family of Hamiltonians (1.1) to the ones with trace one : 



tr 



d 2 Hp 
dA 2 



= 1. 



(1.2) 



The idea is to set up a transformation 1Z that maps 
a Hamiltonian H into a rescaled Hamiltonian 1Z(H), 
such that irrelevant degrees of freedom are eliminated. 
The transformation 1Z should have roughly the following 
properties: TZ has an attractive fixed set (trivial fixed 
set) of integrable Hamiltonians that have a smooth in- 
variant torus with the frequency u>o- Every Hamiltonian 
in its domain of attraction T> has a smooth invariant torus 
with frequency vector u>o- The aim is to show that there 
is another fixed set A which lies on the boundary dT> (the 




critical surface) and that is attractive for every Hamilto- 
nian on &D . 

The transformation 1Z is defined for a fixed frequency 
vector u)o with three incommensurate components. We 
choose u)o — (c 2 ,!?, 1), where a » 1.3247 satisfies a 3 = 
a + 1 (named the spiral mean) . From some of its proper- 
ties, a plays a similar role as the golden mean in the two 
d.o.f. case (34). The analogy comes from the fact that one 
can generate rational approximants by iterating a single 
unimodular matrix N. In what follows, we denote reso- 
nance an element of the sequence {i/k — N x Vi, k > 1} 
where V\ = (1, 0, 0) and 



N 



The word resonance refers to the fact that the small de- 
nominators u>o • Vk appearing in the perturbation series 
or in the KAM iteration, tend to zero geometrically as 
k increases (u>o • Vk = cr 3 ~ k — ► as k — > oo). We no- 
tice that ojo is an eigenvector of N, where N denotes the 
transposed matrix of N. The spectrum of N is composed 
of one real eigenvalue a^ 1 (with ljo as eigenvector), and 
two complex conjugated eigenvalues Ai ± i\ 2 = ^/ae Tia 
(which are of norm larger than one). We denote by 
fi± = fl^ ±ifl( 2 \ the eigenvectors of N associated with 
these complex eigenvalues. 

Our hypothesis (which is also the starting point of a gen- 
eralization of Greene's criterion in Refs. |ll| , p^| ) is that 
the sequence {fk} plays a leading role in the breakup of 
the invariant torus with frequency vector ojq. 
In this paper, we study the extension of the ideas de- 
veloped by Escande and Doveil |35|-|37| to three d.o.f. 
Hamiltonian systems. We build an approximate scheme 
by considering the three main resonances Vi, i/ 2 , and 
1/3. The renormalization focuses on the next smaller 
scale represented by the resonances is 2 , ^3, together with 
1/4 = Nv 3 = 1/1 — 1/3. It includes a partial elimination 
of the perturbation (the part which can be considered 
nonresonant on the smaller scale, namely the mode 
a shift of the resonances, a rescaling of the actions and 
of the energy, and a translation in the action variables. 
It is, in spirit, close to the type of transformations con- 
sidered in Refs. §,§. 

The approximations involved in this scheme are the two 
main ones used by Escande and Doveil: 

a] A quadratic approximation in the actions (as the 
rescaled Hamiltonian 1Z(H) is in general higher than 
quadratic in the actions). 

b] A three resonance approximation: we only keep the 
three main resonances at each iteration of the transfor- 
mation, i.e., we consider the following family of even 
Hamiltonians 

3 

H{A, <p) = H Q {A) + h k (A) cosK • ip), (1.3) 
fe=i 
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where hk denotes the amplitude of the mode of the 
perturbation. 



Hamiltonian (1.1) is isoenergetically nondegenerate if the 
following determinant of order 4 does not vanish 



det 



d 2 H 0H 



dA 2 
dHp 

dA 



dA 




^0. 



On the contrary, if the determinant ( |l.4| ) vanishes, H is 
said to be isoenergetically degenerate. In the following 
section, we define the renormalization transformation for 
both cases. 



II. RENORMALIZATION TRANSFORMATION 

Our transformation is based on the following steps: 
1] We apply a canonical transformation that eliminates 
the first main resonance u%. This is performed by a Lie 
transformation Us ■ (<P,A) l— ¥ if', A'), generated by a 
function S(A, ip). The Hamiltonian expressed in the new 
coordinates is given by 

H> = cxp(S)H = H + {S,H} + {S, H}} + ■ ■ ■ , 

where { , } is the Poisson bracket between two scalar 
functions of the actions and angles: 



{/,<?} 



d£ dg_ 
dip dA 



df_.dg_ 
dA dip' 



and the operator S is defined as SH = {S,H}. De- 
noting by £ the size of hk, the generating function S is 
determined by the requirement that the order 0{e) of the 
mode V\ in H' vanishes: 

{S,H } + hi(A)cos(v 1 -<p) = 0. 

This equation has the solution 

S(A,cp) = S 1 (A)sm{v 1 -<p), 



where 



Si(A) = - 



fti(A) 
w(A) ■ v\ ' 



w(A) 



dH 
dA ' 



This step generates arbitrary orders in the action vari- 
ables [w(A) is linear in the actions]. In order to map 
quadratic Hamiltonians into itself, we expand H' to 
quadratic order in the actions, and we neglect higher or- 
ders. The justification for this approximation is that, as 
the torus is located at A = for H , one can expect that 
for small e, it is close to A = 0. We notice that /12(A) 
and /13(A) are not changed up to order 0(e 3 ). Further- 
more, we neglect all the Fourier modes except 0, 1/2, ^3, 



and 1/4, and all terms of order greater than 3 in e. This 
leads to the expression of H': 

H' = H + h 2 cos(v> 2 • ip) + h 3 cos(v>3 ■ ip) 

+ -({S, hi cos(^i • <p)}) + {S, h 3 cos(iv 3 • ip)}, (2.1) 



(1.4) where () denotes the mean value defined as 



(h)(A)= / h(A,<p) 



T 3 



d 3 y 
(2tt) 3 



The last term of Eq. (2.1) contains the Fourier mode 
1/4 = v\ — 1/3 of amplitude 



MA) -i 



dSj_ 
~dA 



h 3 v 3 ■ -± . (2.2) 



We expand /14 to quadratic order in the actions. 
2] From Eq. (2.1), the mean value term ({S,hiCOs(ui ■ 
ip)}) produces a linear term in the actions. In order that 
the mean value of the linear term in W becomes <jJ ■ A, 
we eliminate this term by a translation in the actions 
A 1 > A + a, where a is of order 0(e 2 ) (so it does not 
produce any other effect up to the second order in e). 
3] We shift the resonances Vk l ~ > Vk-i'- We require that 
the new angles satisfy cos(i/fc+i • ip) = cos(vk • ¥>')> f° r 
k = 1,2,3. This is performed by the linear canonical 
transformation 

(A 1 i P )^(N- 1 A 1 Ni P ). 

We notice that N is an integer matrix with determinant 
one. Therefore, this transformation preserves the T 3 - 
structure of the angles. This step changes the frequency 
uj into Nu)q = (T _1 u;o (since u: is an eigenvector of N 
by construction). 

4] We rescale the energy (or equivalently the time) by a 
factor cr, in order to keep the frequency fixed at u)q. 
5] We rescale the actions: 



H"(A,<p) = \H' [ -, V 



such that Condition (1.2) is satisfied for H". This nor- 



malization condition is essential to the convergence of the 
transformation. 

Similar type of approximate renormalization transforma- 
tions has been defined in Ref. J3^]. The main difference 
is that they used a normalization condition such that the 
Hessian ma trix d 2 Ho/dA 2 is of rank 2, instead of Con- 
dition (1.2). Below, we make the distinction between 
degenerate and nondegenerate Hamiltonians: we explicit 
the renormalization transformation in both cases. 



A. Isoenergetically degenerate Hamiltonians 

The renormalization transformation described in this 
section was derived in Ref. |T^] . The integrable part Ho 
is given by 
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H (A)=u a -A + -(fi-A) 2 1 (2.3) 

where fl is a free vector of norm one: ||f2|| = 

(|^i| 2 + \Q 2 \ 2 + |^3| 2 ) 1/2 = 1. In that case, the Hessian 
matrix d 2 Ho/dA 2 is of rank one (proportional to the 
projection operator on t he H-dircction); thus the isoen- 
ergetic determinant ( [Tj| ) is zero. The relevant direction 
(where there is a twist) in the actions is fi. We expand 
hk(A) in the (ft • ^4)-variable: 

h k (A) = f k +g k fi- A+^m k (ft- A) 2 . (2.4) 



We rewrite the meanvalue terms in H' of Eq. (2.1) as 
Hq(A) + i({5,/ncos(i/i •¥>)}) 

= H (A) + \u x ■ A (Sifci) 

= w • A + aft ■ A + -(1 + /i)(n • A) 2 + const. 

The linear term aft ■ A is eliminated by a translation in 
the actions A' = A + fta/(l + fi). 

The shift of the resonances (Step 3) changes the vector 
fi into Nft. In order to keep a unit norm, we define the 
image of ft by 



ft' 



Nft 

\Nft\ 



(2.5) 



The quadratic term of the integrable part of H' becomes 
cr||7Vn|| 2 (l + fi)(fl' ■ A) 2 j2. We rescale the actions (Step 
5) by a factor 



A = a||iVT2|| 2 (l + M ), 
such that Hq is mapped into 



H^A)=ij a -A+-(ft' -A) 2 



(2.6) 



with ft' given by Eq. ( |2.5| ). The transformation is thus 
equivalent to a mapping acting on an 1 1-dimensional 
space (recall that ft and ft' have unit norm) 

({fk,9k,m k } k=ia x,^) ^ ({/fe)fl , L' m fe}fe=l,2,3;fi') ) 
defined by the following relations 

f' k = a 2 \\Nfi\\ 2 (l + p)f k+1 , 



g' k = a\\Nfl\\g k+1 
1 

-m k+1 



1 + M 

/^=<7 2 ||iVf2|| 2 (l+M)4 0) » 
g' 3 = aWNfiWh^, 

™' - 2 -h {2) 

1 1 A i 



for k = 1,2 



1 + /it 



(2.7) 
(2.8) 

(2.9) 

(2.10) 
(2.11) 

(2.12) 



wher e hf 1 is the coefficient in (ft ■ A) 1 of h± given by Eq. 
\2h. Denoting by 



ft ■ V\ 

UJ ■ Vi 



and 



03 = 



ft ■ Vj, 
UJ ■ Vi 



fix 

7^2 



fh 

a 2 ' 



we obtain explicit expressions for p, and of the renor- 
malization map: 

3 

M = ^x{gi - Pih)(Pigi - Plh - mi), 
4 0) = 



[/3 3 / 3 (5i-/3i/i)+/3i/i33], 

[((A + /3 3 )<?3-2/3 1 /3 3 /3)(0i-/3i/i) 
+(3ifim 3 +(3 3 f 3 m 1 ] , 

[((A + /3 3 /2)m 3 - A (A + 2/3 3 )ff 3 + 3/3 2 /3 3 / 3 ) 
x(.9i -A/i) 

+ (/3x/2 + /5 3 )m lff3 - 3/3 1 /3 3 mi/ 3 /2] . 



Iterating the renormalization map, Eq. (2__5) reduces to 
a rotation; in fact, as uiq is an eigenvector of N with an 
eigenvalue of norm smaller than one, the u>o-direction of 
the vector $7 is contracted. The renormalization trans- 
formation reduces to a 10-dimensional map, where the 
vector ft rotates in the plane (ft^.ft^). It is thus 
parametrized by an angle 9 defined by 

fi = p(ft il) cos9 + ft (2) sin (9). 

If we choose ft 1 - 1 ' such that 

= (cr~ 1/2 cos a, 1, cr 1/2 cos a), 
ft {2) = {<j- 1/2 sina,0,-CT 1/2 sina), 

where Ai ± 1X2 — a x l 2 e^ la are the two complex conju- 
gated eigenvalues of N, the expression of ft becomes 

fi = p(<r~ 1/2 cos(a - 9),cos9, a 1/2 cos(a + 6)), 

where, since = 1, p = [F(9)]^ 1 ^ 2 , and 

F{6) =a- 1 cos 2 (a-6)+cos 2 6 + cTCOs 2 (a + 6). (2.13) 

The parameters (3\ and /3 3 are expressed as functions of 



Pi = «- 5/2 



cos(a — 9) 
[F(6»)]!/2 ' 



_ 3/2 cos(a + g) 
P3 [F{6)Y/ 2 ' 



The expression of the norm || ATfi|| is given by 

1/2 



\Nft\ 



F(a + 9) 
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B. General quadratic Hamiltonians 

For the most general quadratic Hamiltonians, we con- 
sider the family 

H(A, cp) - f(<p) + [u + g(ip)] -A + ^A- [M + m(<p)]A, 

(2.14) 

where M and m are 3x3 symmetric matrices, and g is a 
vector. The matrix M is assumed to be nonzero (its trace 
is equal to one) and g is not parallel to ojq. In Step 3, the 
vector g is renormalized into Ng. Thus the iterations of 
N converge to the plane defined by fi' 1 ' and fl^ 2 \ i.e., 
the Wo-direction of the perturbation is contracted (as the 
modulus of the eigenvalue associated with lj is lower 
than one). In Hq, except the term u)q ■ A which is kept 
fixed by renormalization, in all the other terms of higher 
order in A, the Wo-direction is also contracted. Notice 
that the quadratic term can be written as 



1 



-A ■ MA 



ij = l,2,3 







(«w . i)(nO') . a), 



where fl^ — u>o, and fiW anc ] fti 2 ) are the real and 
imaginary part of the eigenvectors fl± of N. Then it is 
sufficient to consider perturbations which only depend on 
the variables (fiW ■ A) and (fl^ ■ A), and an integrable 
part Hq of the form 



H (A) = uq ■ A 



Aid) 



(«w • A)(n« - a), 



(2.15) 



where Sl^' is a fixed vector (real or imaginary part of the 
complex eigenvectors of N), and the matrix mo, with el- 

(i i) . 

ements m , is symmetric and will be a variable of the 
renormalization map with the restriction that the trace of 
the Hessian matrix is equal to one. There are two direc- 
tions f2W and fi( 2 J of twist in the actions. The Hessian 
matrix is no ninv ertible, but, in general, the isoenergetic 
determinant (1.4) is nonzero. We write h k , for k = 1,2, 3, 
in the (fiw . A)-variables: 

h k (A) = f k +J2 ff£°(n (i) -A) 

t=l,2 

+i ^ ^'(fl^'AJ^.A), (2.16) 



*,J'=1,2 



where the matrices m^, whose elements are mi, , are 



sym metri c. We expand /14 given by Eq. ( |2.2| ) such that 
Eq. ( |2.16] ) defines also the coefficients of the Taylor ex- 
pansion of hi. The shift of the resonances (Step 3) 
changes Vt^ and Q,^ into 

iVO« = A^ 1 ) - A 2 ^ 2 \ 

nq. {2) = x 2 n (1) + Ain (2) . 



This is equivalent to a rotation in the plane (fiW, Jl' 2 ') 
combined with an amplification by a factor (Af-I-A 2 ,) 1 / 2 = 
yfa. This step changes the matrix m k into m' k whose 
elements are 



n'(bi) _ \2_(i,i) 



( X < 2 ) _L \2_(2.2) 



A|m^' i; + 2AiA 2 mjf + X z 2 m y k 



'(1,2) '(2,1) 

m fe = m fe 



= -A 1 A 2 mi 1 ' 1) + (A? - XfimX"' + X 1 X 2 m k ^' , 
TO f 2 ) = \lm^ - 2A 1 A 2 mi 1 ' 2 ) + A 2 mi 2 ' 2 \ 

The matrix mo is changed into m' according to the same 
formulae. The vector g k = {g k ,g k 2 ^) is mapped into 
9k = (9k > 9k 2 " 1 ) wnose elements are 

'(!) _ \ J 1 ) , \ „(2) 

9k - A i3fe + x 2 g k , 



(1,2) 



,(2,2) 



1(2) 

9 k 



\ n W 4- X „ (2) 



The rotation (under the action of N) of g k and mu is 
analogous to the rotation of fi [see Eq. (|2.5| )1; the ampli- 
fication is compensated by the rescaling of the actions to 
avoid divergences of the transformation. We rewrite the 
mean- value term in H' as 

({S, h t cos(i/i • ¥>)}) = const + ^ a W ^ (i) ■ A 

i=l,2 

+ E • A)(ft w • A). 

The linear terms a^fl" • A are eliminated by a trans- 
lation in the actions. The mean-value of the quadratic 
part of H 1 is £^=1,2 <r(m^ j) + ^(^)(n w • A)(ft (j) • 
A)/2, and thus the new Hessian matrix is given by 



E»,j=i,2 cr(m (l j) + ^ & )n w <g> rj (j) , where the elements 

of the matrix n (i) <g> are (fi (l) <g> f2 (i) ) fe/ = O^fy 00 . 
In order to have the trace of the Hessian matrix of the 
rescaled Hamiltonian equal to one, we rescale the actions 
(Step 5) by a factor 



X = a 



/ , v""o 

,3=1,2 



(m (lj) + /i (i ' j) )tr(n (i > <g> n^). 



The approximate transformation is equivalent to a map- 
ping acting on a 20-dimensional space (since the matrices 
m are symmetric and mo has a constant trace): 

{{Jk,9 k i m k ' m }fe=l,2,3;i,j=l,2) 

defined by the following relations 

f'k = ^fk+l, 
nli) lU) 

9 k =(X 9 k +v 

"(iJ) a 

m k = jmk+i, 
u(i,j) a lUJ) 

m ° = A m ° 
for k = 1, 2, 3 and i, j = 1, 2. 
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III. RENORMALIZATION FLOW 



For each scheme (Sees. II A and II B), the numerical 
implementation shows that there are two main domains 
separated by a critical surface: one where the iteration 
converges towards a family of integrable Hamiltonians 
(trivial fixed set), and the other where it diverges to in- 
finity. 



A. Reduction to degenerate Hamiltonians 

As a first result, we numerically observe that the trans 
formation a cting on nondegenerate Hamiltonians consid- 



ered in Sec. [I B tends to the degenerate ones of Sec. [I A , 
on the critical surf ace. More precisely, if we consider 
Hamiltonians (2.14) with M of rank 3, th e co ntraction 



in the u;o-direction, as explained in Sec. [IB, reduces 



the rank of M by one; the 3x3 matrix M is thus re- 
duced to a 2 x 2 matrix too. Furthermore the numerical 
results show that the renormalization reduces this rank 
to one when we iterate on the critical surface. Figure 1 
shows the evolution, under the renormalization map, of 
the determinant of to . The upper curve corresponds to 
a starting Hamiltonian in the domain of attraction of the 
trivial fixed set, and the lower one corresponds to itera- 
tions on the critical surface (both evolutions start with 
the same quadratic part). We also check that the deter- 
minant of the matrices rrik tends to zero. Furthermore, 
the directions of the vectorial parameters (g a nd to ) tend 
to be aligned by t he i teration : Hamiltoni ans ( 2.16|) tend 
to Hamiltonia ns (2.4), and Hamiltonian (2.15) tends to 
Hamiltonian (2.3), all with a same direction fi. In or- 
der to characterize this, we define fi as the unit vector 
with the same direction as g\. We compute the norm of 
r?2 — and r?3 — fl, where f?2 (resp. O3) is a unit vector 
with direction gi (resp. 93). Moreover, in order to see 
that the quadratic terms are proportional to ft ® fi, we 
compute the norm of too — cl~i <g) $1 where c is defined by 
the constant trace of too (similar calculations have been 
done for the other matrices TOfc). These differences tend 
to zero as we iterate a Hamiltonian on the critical surface 
(they are of order 10~ 5 after 20 iterations on the critical 
surface). 

From these observations, we conclude that isoenergeti- 
cally degenerate and nondegenerate Hamiltonians belong 
to the same universality class. According to the general 
renormalization-group picture, the corresponding critical 
invariant tori are predicted to have the same type of scal- 
ing properties. 

We lack an explanation of the mechanism of this second 
reduction of the rank. We remark that this second re- 
duction is not just an effect of the rescaling (steps 3 to 
5) as it is the case for the u>o-contraction : the second 
reduction does not happen outside the critical surface. 
We conjecture that this reduction will also occur in an 
exact renormalization scheme, but this point has not yet 



been explored. 

We remark that the choice of the normalization condi- 
tion (1.2) seems essential to obtain a nontrivial attrac- 
tor. Other choices [B2|, like det mo = 1, do not lead to a 
critical attractor and the iterations appear to diverge on 
the critical surface (one eigenvalue of too tends to zero 
and the other one to infinity). 

The reduction to rank one allows us to work, for the pre- 
cise analysis of the attrac tor, with the data obtained for 
the degenerate case (Sec. HA). 



B. Trivial attractor 

The domain of attraction of the trivial fixed set is the 
domain where the perturbation of the iterated Hamilto- 
nians tends to zero. However, the renormalization tra- 
jectories in this domain do not converge to a fixed Hamil- 
tonian but converge to a smooth quasiperiodic set of in- 
tegrable Ham ilton ians. This can be explained by looking 
at the map (2.5). The eigenvalues of N are a~ 1 and 
v^e ±M where a « 2tt x 0.3880 (a = arccos(-cr 3 / 2 /2)). 
The map (2.5) leads asymptotically to a rotation of an- 
gle a in the plane (nW, J^ 2 -*), afte r a co ntraction in the 
u;Q-direction, as e xplained in Sec. II A. The values of 



the rescalings (2.6) at the trivial fixed set are given by a 
smooth function of 9. It is given explicitely by 



A(0) = a 



2 F(a 



F(9) 



where F is given by Eq. ( |2.13 ). This trivial rescaling 
curve is depicted in Fig. 2. Since a/2n is close to 7/18, 
the evolution of A oscillates approximately with period 
18. 



C. Critical attractor 

On the critical surface, the renormalization flow con- 
verges to an attracting set. This set has a codimension 
1 stable manifold, i.e., one expansive direction trans- 
verse to the critical surface. This set plays, for the sys- 
tem we consider, the same role as the nontrivial fixed 
point of the renormalization-group transformation for 
quadratic irrational frequencies in two d.o.f. Hamiltonian 
systems. In particular, its existence implies universality 
for one-parameter families crossing the critical surface. 
Different trajectories of the transformation display the 
same values of the rescalings with a different order but 
with a universal statistical distribution. We define expo- 
nents that characterize the universality class associated 
with the spiral mean. The mean-rescaling is defined by 

A = linin^oo ( A; ) > wnere <\? is the value of the 
rescaling after j iterations on the critical surface. We 
also calculate the largest Lyapunov exponent k which 
measures the approach to criticality as a function of the 



G 



coupling constant, of the universality class. The result we 
found is that these limits do not depend on the Hamilto- 
nian on the critical surface where we start the iteration 
nor on the initial choice of i~2. The coefficients k and A 
depend only on ujq. Numerically, we find k ~ 0.6427 and 
A w 3.1479. 

We provide numerical evidence that this attract or is 
s tran ge and nonchaotic. We remark that Eqs. (2.7)- 
(2.12) of the renormalization map have the form of a 
nonlinear system {//., gk-,mk} (for k = 1,2,3) driven by 
a quasipcriodic variable 9 t— > 9 + a, give n by the evolution 



of the vector fl represented by Eq. (2.5). This type of sys- 
tems has been analyzed in Refs. [JL8- 



22 , where one of the 



main conclusion is the existence of strange nonchaotic at- 
tractors. In order to analyze the nontrivial attractor from 
this perspective, we show in Fig. 2 a two-dimensional 
plot of the time series (9j,Xj) of the scaling factor Xj, 
and the angle 6j, where j is the index of the iteration 
on the attractor. This figure shows that A appears to 
be a continuous (one to one) function of 9. The evolu- 
tion of the rescaling Xj displays an approximate period 
18 behavior, similar to the one observed on the trivial 
attractor. We remark that for the trivial attractor X(9) 
is smooth, while for the critical attractor X(9) has a set of 
cusps (nondifferentiable points). Since the driving map 
9j+i = 9j + a fills the circle densely, and the renormaliza- 
tion map is smooth, the function X(9) must have a dense 
set of cusps. In Fig. 3, we show the corresponding plot 
for the parameter g\ . This picture is the renormalization 
trajectory of a single initial Hamiltonian on the attrac- 
tor. It shows that gi(9) is not a single valued function of 
9. Similar pictures are obtained for the other coordinates 
of the map. 

In order to analyze the structure of the attractor in more 
detail, we consider it from a different point of view: We 
take a set of initial conditions on the critical surface (not 
on the attractor) parametrized by an angle 9 varying over 
a small interval [9i, 62]. Figure 4a shows the projection of 
this segment on the plane (9, g%). Figure 4b shows the im- 
age of the projection after 100 iterations, and Fig. 4c af- 
ter 350 iterations. This shows that as the segment comes 
closer to the attractor the number of steep oscillations 
becomes larger, suggesting that in the limit they corre- 
spond to discontinuities of the attractor. This behavior 
is similar to the observations of Ref. Il8] and is compat- 
ible with the results of Ref. |^|. These oscillations are 
associated to abrupt changes of the signs of the coordi- 
nates. One can conjecture that on the attractor, there 
is an infinite number of such changes of sign. The same 
kind of phenomena is observed for the other coordinates. 
In order to see the effect of these changes of sign in the 
renormalization dynamics on the attractor, we compare 
in Fig. 5 the power spectrum of the time series of gi, with 
the one of \gi\. The second one indicates that the behav- 
ior of \gi\ is very close to quasiperiodic (with frequencies 
1 and a/27r), while g\ shows a broad spectrum. 
A further insight into the structure of the attractor can 
be obtained by following the trajectories for a single ini- 



tial 9 and a set of different choices of the other coor- 
dinates at a given "time" (i.e., after a fixed number of 
iterations of the map). After the relaxation transient, 
one observes that all initial points fall into one of eight 
points of the attractor (in the example of Ref. 1 1^] , there 
are two such limiting points). The iteration of these eight 
points generates eight branches of the attractor that are 
not completely disjoint. In order to visualize these eight 
branches, we display in Fig. 6 the values of an observable 
B that allows to distinguish them. A suitable choice is 
the weighted sum of the coordinates of x — {x\, . . . , Xg) 



with 



i=l 



N 



where Xi(j) is the j-th iterate of an initial condition of the 
coordinate Xi. Thes e bran ches are related by symmetry, 
as described in Sec. 



hid 



In summary, the numerical results suggest the following 
characterization of the attractor: It is a set composed of 
eight branches that differ only in the signs of one or more 
hk (the amplitude of the Fourier mode u^)- Each of the 
branches has a dense set of discontinuities. It seems to 
be an example in 9 dimensions of the type of attractors 
described in Ref. fTq| . 

In order to characterize the singularities, we com- 
pute the correlation dimension of the attractor accord- 
ing to the method developed by Grassberger and Pro- 
caccia [[58 39 . We did not find a clear result, but there is 
some evidence that the dimension has a value around 1, 
in agreement with Figs. 3 and 6, and with the projection 
of the attractor on the plane (g±, mi) depicted in Fig. 7. 



D. Symmetries of the transformation 

In this paper, we have considered only even perturba- 
tions, i.e., such that the Fourier modes are only cosine 
terms. The results can be extended to the more general 
case, including non-even modes v>k (k — 1,2,3), by the 
following symmetry arguments (4fjj^]. These arguments 
allow also to understand the relation between the eight 
branches of the attractor. 

We analyze the effect of a shift of the origin of the angles 
on the renormalization transformation. We denote this 
shift as 



T, 



: (fi h-> ip 



6. 



The KAM transformation (Step 1) commutes with Tg. 
The action of Tg on the shift of the resonances (Step 3) is 
characterized by the following intertwinning relation J4l| : 



7 



TZoT g =T^ e oTZ, 



(3.1) 



where 1Z denotes the renormalization transformation. 
Applying this relation to the critical attractor gives 
the relation between the renormalization trajectories for 
Hamiltonians H(A, (p) and the ones for Hamiltonians 
H(A, ip + 6). With this type of shift, any Hamiltonian 
containing only the three modes f^, k = 1,2,3 (with sine 
and cosine terms), can be put into a cosine representa- 
tion. Thus the attractors for these models are directly 
linked by symmetries to the attractor found in the even 
case. 

The eight branches of the critical attractor are mapped 
into each other by symmetries of this type, realized by 
shifts in the origin of the angles by 6k = nu^. This cor- 
responds to the eight possible choices of the signs of the 
three modes. 



IV. CONCLUSION 

This paper provides numerical results indicating that 
for the spiral mean frequency vector, the critical surface 
of the approximate renormalization transformation is the 
codimcnsion 1 stable manifold of a strange nonchaotic 
attractor. This feature depends strongly on the charac- 
teristics of the eigenvalues of N. Moreover, the numerical 
results suggest that the renormalization transformation 
can be reduced to an isoenergetically degenerate fam- 
ily of Hamiltonians at criticality. These remarks give 
new insights for the setup of a systematic renormaliza- 
tion transformation, in the spirit of Refs. pd|,p|,|9| . 
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FIG. 1. Evolution of the determinant D„ of the matrix mo 
(as a function of the number of iterations n) : the continuous 
line is for a trajectory of the renormalization transformation 
on the critical surface, and the dashed line is for a trajectory 
inside the domain of attraction of the trivial fixed set. 



10 



X(9) 







i ■¥ E 

i ; i • ?! 

7 toil 




\ M 

W v. 


■ w 



-TT 6 



FIG. 2. Values of the rescalings A as a function of the angle 
6 between ft and f2 (1) in the (fi (1) , fi (2) )-plane. The regular 
curve is for the trivial fixed set, and the singular curve is for 
the nontrivial fixed set. 
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FIG. 3. Values of g\ as a function of the angle 9, on the 
critical attractor of the renormalization map. 
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FIG. 5. Power spectrum of the evolution on the nontrivial 
fixed set, (a) of the coordinate gi, and (b) of the absolute 
values \gi\. 
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FIG. 4. Evolution of a set of initial conditions on the criti- 
cal surface, parametrized by an angle 9 varying over an inter- 
val ^2]. We depict the projection of this set on the plane 
(9,gi) : (a) initial set, (b) after 100 iterations, and (c) after 
350 iterations. 
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FIG. 6. Values of an observable B as a function of 9, on the 
critical attractor: it shows the eight branches of the attractor. 
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FIG. 7. Projection on the plane (gi,mi) of the critical at- 
tractor of the renormalization map. 
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